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We discuss a proposal for the efficient generation of the 
maximally entangled atomic N-GHZ state in a spinor-1 con- 
densate by driving internal state atomic Raman transitions 
using (classical) laser fields. We illustrate the dynamics in 
terms of a quantum Zeno subspace, and identify the resultant 
atomic elastic collision in facilitating the deterministic entan- 
glement creation. Our proposal can be readily implemented in 
several laboratories where ferromagnetic spinor condensates 
(of 87 Rb atoms) are investigated. 
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The continued success of the research on quantum de- 
generate atomic gases has led to intense exploration of 
their potential application for quantum information sci- 
ence. A substantial topic of recent interest is the creation 
of spin squeezed atomic states that display multi-particle 
inseparable correlations [1,2]. While it is debatable how 
such inseparable quantum correlations among identical 
particles are related to entanglement between distinct 
parties [3,4], it has nevertheless been fruitful to inves- 
tigate these unexplored territories. Maximally entangled 
states enjoy special attention; in addition to being states 
of maximum inseparable correlations or entanglement, 
they compose of a special class with which theoretical 
discussions can become most transparent. 

Recently, a deterministic protocol was suggested [5] 
for creating maximally entangled pairs, triplets, quar- 
tiles, and other clusters of Bose condensed atoms starting 
from a condensate in the Mott insulator state [6]. The 
work of [5] concerns a single optical well and involved 
a small number of atoms. In this limit, the evolution 
dynamics of the system is analytically solvable, which al- 
lows the appropriate timing of the external laser fields 
(in comparison with atom elastic collision strength) to 
be identified for creating maximal atomic entanglement. 

This article clarifies the operating mechanism of the 
seemingly simple protocol [5]. More importantly, we in- 
vestigate its application to large numbers of condensed 
atoms. To our surprise, the few atom protocol with 
classical laser fields driving atomic Raman transitions 
remains efficient in generating massive atomic entangle- 
ment of condensed atoms. The answer, it turns out, can 
be obtained rather conveniently understood in terms of 
a quantum Zeno subspace [7] . Within this subspace, the 
effective atom interaction and the Raman dynamics re- 
duces to a well-known form studied earlier [8]. What is 
crucial for creating maximally entangled atomic states is 



to start from a condensate with a fixed number of atoms, 
(be it even or odd [8]). 

This paper is organized as follows. First, we briefly in- 
troduce our system and the setup for generating maximal 
entangled atomic states. Then we provide a step by step 
illustration of the underlying mechanism for the creation 
of such massively correlated atomic states. We briefly re- 
view the concept of a quantum Zeno subspace [7] . This 
is followed by a review of a spin-1 condensate, paying 
special attention to the case of ferromagnetic atom-atom 
interactions (as for 87 Rb atoms). We then discuss an in- 
teresting decomposition of the interaction Hamiltonian 
into three sub-SU(2) spaces as first suggested in [9], Fi- 
nally, we show that the external laser fields in the ap- 
propriate Raman configuration create a quantum Zeno 
subspace, within which the effective interaction is of the 
desired form for generating spin squeezing and maximal 
atomic entanglement. We conclude with some illustrat- 
ing numerical results. 

Our main result can be summarized in Fig. 1 below. 
We consider a spin-1 [F = 1 atoms with three Zeeman 
sub-levels) condensate with a fixed number (N) of atoms 
that are all in the |+) = \M F = 1) or |-) = \M F = -1) 
state initially. We claim that if a strong off-resonant 
Raman transition is established between the |+) and |— ) 
state, then at predictable instants of the time evolution, 
the condensate becomes maximally entangled, i.e. its 
atomic internal state becomes an N-atom GHZ state oc 
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FIG. 1. Illustration of our basic result. 

Quantum Zeno subspace The quantum measure- 
ment process itself can strongly affect the dynamics of 
a system. An often cited example in this regard is the 
famous quantum Zeno effect [10,11]: a watched pot does 
not boil, or in a physical system, frequent measurements 
of the initial state hinders its dynamics to evolve to other 
states. Recently, Facchi and Pascazio provided new in- 
sights into the Zeno effect, by reformulating it in terms 



of an adiabatic theorem within different sectors of the so- 
called quantum Zeno subspaces [7] , which are closely re- 
lated to the powerful decoherence free subspace in quan- 
tum computing studies [12,13]. In terms of the notation 
of [7], the concept of a quantum Zeno subspace can be 
illustrated by a system Hamiltonian 

H n =H + SlH mcas . (1) 

In the limit of an infinitely large parameter a, the time 
evolution of Ha is diagonal with respect to spaces of 
Off mcas , i.e. for any Hamiltonian of the form (1), the 
largest term aH mcas decomposes the complete Hilbert 
space into different subspaces or quantum Zeno sub- 
spaces. The dynamics associated with the small term 
H is slaved to slow adiabatic evolutions within each of 
these orthogonal subspace. 

Spinor condensate For ferromagnetic interactions 
(e.g. 87 Rb), as was shown before [14], the spatial mode 
of different condensate component 4>j{f) (normalized to 
unity) is identically the same <f>(r). Thus, a spin-1 con- 
densate can be described with a J=+ , ,- ( a j=+,o.-); tnc 
single atom annihilation (creation) operator for each spin 
component state. Apart from constant terms that de- 
pend on the total number of atoms N = N + + N + iV_ , 
(Nj = Ojdj), the atom- atom interaction becomes [15-17] 

H = uL 2 , (2) 

where the coefficient is u — c 2 J \<f)(r)\ 4 dr/2 (u < for 
ferromagnetic interactions). c 2 = 4irh 2 (a2 — a )/3M with 
a2 and ao the respective scattering lengths in the to- 
tal spin F = 2 and channels of two colliding spin-1 
atoms. The pseudo angular momentum operator L is 
defined according to the Schwinger representation, and 
with its components given by 

L + = L x + iL y = V2(a} + av + aja_), 

L_= , L z = a) + a + — a^a_. (3) 

It is well known [9] that such an effective interaction 
(2) is formally SU(2) symmetric, thus would have not 
generated any nonlinear dynamics if the L/s were the 
SU(2) Schwinger representation of two bosonic modes 
[18]. For the spin-1 condensate being considered here, as 
was emphasized before [9], the pseudo angular momen- 
tum operator L does NOT satisfy the Casimir relation 
[L 2 7^ N(N + 1)], and therefore can not really be con- 
sidered an angular momentum operator [9]. Using the 
Gell-Mann decomposition of SU(3) into three SU(2) sub- 
spaces U, V, and T, it was found [9] that 

L 2 = AT 2 + 2(V+U+ + V-U-) 

+ (N-e + )(N-e-)/2-2(Y-Y ) 2 , (4) 

with e± = -3/2 ± V2 and Y = -N/6 - 1/4. It is 
important to note that operators belonging to different 



SU(2) subspaces do not always commute. We adopt the 
convention that 

T+ = a\a-, T z = ^{a\a+ - ala_), (5) 

( U + + ) = a ± a °' ( u* z ) = \( a ± a± ~ a o a °)' ( 6 ) 

and the hypercharge Y = (N + + iV_ — 2A^o)/3, which 
physically corresponds to the quadrupolc moment of 
atomic populations Nj— +t o,-- 

Effective interaction and maximal entangle- 
ment generation A Hamiltonian of the form T 2 [the 
first term in Eq. (4)] can be used to generate maximally 
entangled or spin squeezed states [1,2,8,19]. To selec- 
tively suppress other terms in the Hamiltonian Eq. (2), 
we resort to the idea of quantum Zeno subspace as re- 
viewed earlier. In this case, to limit the dynamics to 
the manifold composed of internal states |+) and |— ), 
we simply Raman couple them with strong external laser 
fields, in the form f2(T + — T_)/2i = QT y as shown in Fig. 
1. In the spirit of the adiabatic theorem within the Zeno 
subspace [7] , our system dynamics is now governed by 

h' = 4uT 2 + nr y , (7) 

an effective Hamiltonian that is known to generate max- 
imally entangled states [5,8]. 

The shortest time for creating a maximally entangled 
state when initially all atoms are in either |+) or |— ) 
is, 7r/4|u|, a factor of two longer than that from a bare 
interaction of the form H' = 4uT 2 [5]. This point be- 
comes clear when we transfer to the interaction picture 
\\p{t))i = exp(iQT y t)\ip(t)) as now SlT y becomes the 
largest term in the system Hamiltonian. Equation (7) 
arises when the limit il ^> N\u\ is satisfied. The effective 
Hamiltonian in the interaction picture then becomes 

H! nt = e inT »*(4uT z 2 ) e - inT «'* 

= Au(T z cos Clt - T x sin at) 2 

w 2u(T 2 + T 2 ) = 2uT 2 - 2uT 2 , (8) 

where the approximation is due to averaging over the 
rapid Rabi oscillations. The first term oc T 2 in Eq. (8) is 
SU(2) symmetric in the two dimensional subspace of |+) 
and |— ), thus it can be neglected as it only introduces 
an overall phase factor. The second term, cx T 2 (with 
a coefficient a factor of 2 smaller), is similar to the T 2 
term [5,8], thus explaining the dynamic creation of max- 
imally entangled states. The entanglement remains for 
the Schrodinger picture based atomic states as the two 
pictures coincide at integer periods of the Rabi oscilla- 
tion. The rapid time average results in a factor of two 
reduction of the effective interaction strength. 

We now present selective results from numerical simu- 
lations based on the complete Hamiltonian (2) plus the 
Raman term aT y and compare with the approximate re- 
sults based on the adiabatic Hamiltonian (7). 



Figure 2 demonstrates the validity of a quantum Zeno 
subspace (|+) and |— )) with increasing values of the effec- 
tive Rabi frequency and for two different atom numbers. 
Indeed, when becomes larger as compared to N\u\, 
the populations become increasingly localized within the 
isospin subspace (|+), |— )), and eventually the Hamilto- 
nian H' correctly describes the dynamics and maximally 
entangled N-GHZ states are created at t = n/4\u\. The 
creation of a Zeno subspace can be further illustrated by 
watching the time dependent population distributions as 
in Figs. 3 and 4. 
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FIG. 2. The average fraction of atoms in the Zeno subspace 
n z = (N+ + N-)/N (solid line) and in state jO), n = {N )/N 
(dashed line) for N — 10 (left panel) and 50 (right panel) at 
t = 7r/4|u| = 5 (ms). The qualitative dependence is similar 
at other times. u/(2tv) — —25 (Hz) and initially all atoms in 
state +). 

We found that if a strong Raman coupling is estab- 
lished between states |+) (|— )) and |0), a quantum Zeno 
subspace composed of atoms in these two states is also 
created. However, maximally entangled states do not oc- 
cur as the projection of L 2 into this space does not yield 
the required nonlinear interaction V 2 (U 2 ) [9]. 




FIG. 3. The time dependence of (N+) for TV = 20, 
u/(2n) = —25 (Hz), and initially all atoms in state |+). The 
solid line is from the complete Hamiltonian while the dashed 
line is from the adiabatic Hamiltonian in the Zeno subspace. 
At n = 40N\u\ = 40 x (20|u|) = (2tt)20 (kHz) (lower panel), 
the Zeno subspace is well established. 



We now discuss the generation of N-GHZ states when 
f2 S> iV|u|. Because the interaction picture Hamiltonian 
(8) is similar to that of Molmer's model as in [8], we 
anticipate similar dynamical behaviors. We adopt the 
notion of the maximal entangled fraction and define the 
optimized overlap 



F max = ma,x(\ N (E s ^^\ijj)\ 2 ) 

e,(j>,n 

= 0.5(|a| + |/3|) 2 , 



(9) 



of the N-atom wave function with an entangled 

coherent spin state [20] 
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(10) 



with0 = 7T-0, <j> = ir + (j>, a = (9 } (f>\ip), and (3 = (0,<j>\ip). 
For F max we note rj = arg(/3) - arg(a). 
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FIG. 4. The same as in Fig. 3 but the probability for 
all atoms in |+). Cjv = {ip(t)\a^ \vac) /VW.. Note the the 
creation of a N-GHZ state at time ~ 5 (ms) (lower panel) . 




FIG. 5. The time dependence of F max (solid line), 
\a\ 2 (dashed line), and j/3j 2 (dotted line) for N = 20, 
m/(2tt) = -25 (Hz), and fi = (2tt)10 (kHz). 



In Fig. 5, we display the time dependence of -F max with 
the use of the adiabatic Hamiltonian (7). We see two 
interesting regimes. For times up to l/(N\u\), the wave 
function is almost a coherent spin state, so the mean field 
theory can be used to describe the system dynamics [19]. 
Around t — 7r/4|u| = 5 (ms), however, we obtain a high 
fidelity entangled coherent spin state, where the mean 
field theory completely breaks down. We also find that 
the angles # max , max , and ?7 max are essentially constant 
in the immediate neighborhood of t = n/4\u\ = 5 (ms) 
(in the interaction picture), thus a phase sensitive detec- 
tion of the massive entanglement can be accomplished 
without extremely precise timing [19,21]. F max has been 
proven to be a useful measure of N-atom entanglement: a 
N-atom pure state \ip) is entangled if F max > 1/2 (above 
the dash dotted line in Fig. 5) [22,23]. 

Finally, we note that the maximally entangled state 
1+)®^ + |— )® JV as generated with our method is sta- 
ble against elastic two body collisions, thus is ideal for 
studying multiple atom entanglement. Furthermore, our 
protocol works, even when there is a non-zero magnetic 
field that breaks the level degeneracy for different spin 
components. The Zeeman interaction gives rise to 

Hb = —TiujL(a) + a + — ala_) + M(a^ + a + + a^a_), (11) 

where the linear Zeeman term is proportional to the Lar- 
mor precessing frequency ojl = BfiB/h, with [1b the 
magnetic dipole moment for state |±), and the quadratic 
term is proportional to S cx B 2 , one half the energy re- 
leased from the generation of two atoms in state |0) due 
to the elastic collision of one atom in state |+) with an- 
other one in state |— ). In the quantum Zcno subspace, 
a^a + + a)_ci- = N, remains a constant, so the quadratic 
Zeeman term does not affect the dynamics of generating 
maximally entangled N-GHZ states. The linear Zeeman 
term is oc T z , which simply induces a twisted phase on 
the state |±). In the rotating frame, the atom-atom in- 
teraction term remains the same 4uT 2 , while the Raman 
coupling £lT y now contains time dependent phase fac- 
tors e ±luJLt . Nevertheless the complete system dynamics 
remain the same. Therefore we expect our theoretical 
protocol to remain effective. 

In conclusion, we have proposed a protocol for robustly 
generating maximally entangled atomic states in a con- 
densate. We have illuminated the operating principle in 
terms of a quantum Zeno subspace. Our protocol can 
be implemented in currently available condensate sys- 
tems with ferromagnetic interactions, where the spatial 
mode functions are identical for each of the components 
[14]. It can also be applied to spin-1 condensates with 
anti- ferromagnetic interactions (e.g. 23 Na atoms [24]), 
which may possess larger values of exchange interaction 
u, thus leading to faster generation of entangled conden- 
sates. Although the different equilibrium mode functions 
in the latter case can lead to modulational instability, 
thus increased decoherence of the maximally entangled 
state [14]. Our protocol starts with all condensed atoms 



in either state |+) or |— ). Provided the time scale l/|u| is 
reasonably short as compared to the inverse of a typical 
condensate phonon frequency, the dynamics for gener- 
ating entangled condensate is also mechanically stable, 
irrespective of whether the |+) or |— ) condensate com- 
ponent is miscible (for anti-ferromagnetic interactions) or 
not (ferromagnetic interactions). 
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